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1 Abstract. An XMM-Newton observation of the cool (kT = 2.1 keV) cluster A1983, at z = 0.044, is presented. 

Gas density and temperature profiles are calculated over the radial range up to 500 h^, 1 kpc, corresponding to 
~ 0.35 V2oo- The outer regions of the surface brightness profile are well described with a j3- model with f3 = 0.74, 
but the central regions require the introduction of a second component. The temperature profile is flat at the 
exterior with a slight dip towards the centre. The total mass profile, calculated from the temperature and density 
information assuming hydrostatic equilibrium, is consistent with an NFW profile, but with a low concentration 
parameter c = 3.75 ± 0.74, which may be due to the cluster not being totally relaxed. Published optical data 
are used to calculate the M/Lb ratio profile and the overall iron mass over luminosity ratio. The M/Lb ratio 
£ — ' profile shows that, at large scale, light traces mass to a reasonable extent, and the M/Lb ratio at 0.35r2oo 

(M/Lb = 135 ± 45 /15O-M0/L©) is consistent with the trends with mass observed in the optical. The iron mass 
over luminosity ratio is about two times less than that observed for a cluster at 5 keV. The gas mass fraction 
rises rapidly in the central regions to level off quickly at ~ 200 h^ kpc; the value at 0.35 r2oo is ~ 8%. The 
scaling properties of the emission measure profile are consistent with the empirical relation M gas oc y 1 - 94 ; use of 
the standard self-similar relation M gas oc T 1,5 results in a scaled profile that is a factor of about two too low as 
compared to the reference mean profile for hot clusters. Comparison of the entropy profile of this cool cluster 
with that of the hot cluster A1413 shows that the two profiles are extremely well scaled using the empirically 
determined relation S oc y - 65 , suggesting that the slope of the S-T relation is shallower than expected in the 
standard self-similar model. The form of the two entropy profiles is remarkably similar, and there is no sign of a 
larger isentropic core in the cooler cluster. These data provide powerful agruments against preheating models. In 
turn, there is now increasing observational support for a trend of / gas with system mass, which may go some way 
, towards explaining the observed scaling behaviour. 
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1. Introduction of a cluster roughly corresponds to S = 200 (e.g., Evrard, 

Metzler & Navarro I1996|l . Assuming a constant gas mass 

In the simp lest m odels of gravitational collapse (e.g., f ract i n / gas = M gas /M v , the total gas mass then scales as 

Bertschingerll985J, supported by numerical simulations Mgag0C T 3/2_ Furt hermore, if the X-ray emission is domi- 

without non-gr avitat ional gas processes (e.g., Navarro, na ^ ed by Bremsstrahlung, then L x ocT 2 . 
Frenk & White fiW7l [NFW] ) and analytical models (e.g., 
Cavaliere, Menci & Tozzi HUMjl . the intrinsic properties 

of galaxy clusters (e.g., X-ray luminosity and tempera- However, it has been known for more than a decade 

ture L X ,T, virial mass My) follow self-similar scaling de- that real clusters do not follow these laws (e.g., Edge & 

rived, assuming HE, from basic virial relations. Such self- Stewart HMD ■ The most notable example is perhaps the 

similarity applies equally to the hot intracluster medium £ X -T relation, with observations suggesting L x oc 7 1 ~ 2 - 9 

(ICM) and the dark matter components. In particular, ( e . g-) Arnaud & Evrard HMfy . In the study of why ob- 

M s oc T 3 / 2 , where M s is the total mass in a sphere of ser ved cluster properties deviate from the expected scal- 

radius R s corresponding to the overall density contrast S, i ngS; it i s the L x -T relation which has been in the spot- 

and T is the cluster X-ray temperature. The virialised part light, mainly because this relation deals with easily observ- 

able global properties. Other important scaling relations 

Send offprint requests to: G. W. Pratt, e-mail: exist, however. Of particular interest for the present work 

gwp@discovery.saclay.cea.fr are the M gas -T and entropy S-T relations, because both 
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of these relations probe variations in gas content and can 
directly be derived from the data. 

Studies of large samples, such as those by Neumann 
& Arnaud (|2001[1 and Mohr, Mathiesen & Evrard l|1999[) 
have shown that, in fact, M gas oc T 19 , a steepening which 
can be interpreted as a dependence of / gas on tempera- 
ture. However, the direct calculation of / gas is fraught with 
difficulties, requiring a robust method for calculating the 
total mass My, and so whether / gas does depend on tem- 
perature is still under debate (e.g., Arnaud & Evrard[l999 
Mohr et al. IT5§9l Roussel, Sadat & Blanchard lMjHI) . 

The entropy is a formidable tool for studying the ther- 
modynamic history of the ICM (e.g., Voit et al 2003). 
Defining the 'entropy' as S — T/nV Z (e.g., Lloyd-Davies, 
Ponman & Cannon 2000), self-similarity implies a simple 
scaling of the entropy with temperature, such that S ex T. 
Ponman, Cannon & Navarro {X22H1) and Lloyd-Davies et 
al. (2000) used ROSAT observations directly to detect en- 
tropies exceeding levels attainable through gravitational 
collapse alone, in the inner regions (0.1 r2oo) of low mass 
systems. From their sample of 20 systems, it appeared that 
the entropy followed the expected self-similar scaling down 
to a certain floor level of 140 hgQ 1 ^ 3 keV cm 2 , implying that 
non-gravitational processes act to set a lower limit to the 
entropy possible for the gas in collapsed haloes. However, 
very recent results with a much larger sample of 66 sys- 
tems (Ponman, Sanderson & Finoguenov 2003), indicate 
that, in fact, there exists a slope in S(T) which is signifi- 
cantly shallower than the self-similar relation throughout, 
such that S oc T~ 65 . This again calls into question the 
fundamental assumption of a constant / gas - 

The commonest explanation for the deviation from 
simple self-similar scaling requires non-gravitational en- 
ergy input, raising the gas to a higher entropy level at 
some time in its history. However, there is little consensus 
regarding the source of this extra entropy. Pre-heating has 
been suggested (originally by Kaiser 119911 and Evrard & 
Henry 1991), although there is no agreement on the as- 
trophysical source responsible for it: early galactic winds 
(e.g., Loewenstein 2000), AGN (Valageas & Silk QUI, 
or possibly both. Other effects may also play a role, like 
radiative cooling (e.g., Muanwong et al. |J001), or vari- 
ation of galaxy formation efficiency with system mass 
(Bryan QUI)- 

Cool clusters are ideal targets for investigating the 
entropy because in such systems, the excess entropy be- 
comes more evident as the processes which produce it af- 
fect the ICM in direct competition with standard grav- 
itational heating. A1983 is well studied in optical, with 
~ 125 galaxy velocities measured in the field by Dressier 
& Shectman (|1988f> . Girardi et al. (1993) derive a velocity 
dispersion of 551+4 7 km s _1 from these data, a relatively 
low value which implies that A1983 is a cool/poor cluster. 

In contrast, A1983 is little studied in X-ray. It was 
observed with Einstein (Jones & Forman 1984), giving 
L x [0.5 - 3.0] keV = 2.12 ± 0.44 x 10 43 erg s" 1 , and has 
been detected in the ROSAT All Sky Survey (Ebeling et 



al. ITggg|> with a quoted L x [0.1 - 2.4 keV] = 0.47 x 
10 44 erg s _1 . These low Lx values again point to the poor 
nature of the cluster. 

This paper first presents results from the new XMM- 
Newton observation of A1983. The entropy and emission 
measure profiles of this cool cluster are then compared 
with those of hot clusters. These are two emphases: a di- 
rect comparison in form between the radial profiles, which 
shows the unparalleled quality of these XMM-Newton 
data; and how the scaling relations involving these quan- 
tities varies across the temperature range between the 
clusters. This paper will show the remarkable similarity 
in form between the profiles which is difficult to recon- 
cile with conventional preheating models, and will use the 
M gas -T and S-T relations to argue that the gas mass 
fraction / gas cannot be constant with cluster mass. 

Taking Ho = 50 km s _1 Mpc~ 4 and 00 = 1, 1' corre- 
sponds to 71.4 kpc at the cluster redshift of 0.044. 

2. Data preparation and analysis 

2.1. Preliminaries 

The XMM-Newton observations of A1983 were taken in 
Revolution 400, on 2002 February 14, and have a total 
exposure time of ~ 32 ks. Pipeline products provided by 
the XMM-Newton SOC, consisting of calibrated event files 
preprocessed with the SAS, are used in this analysis. All 
camera exposures used the THIN1 filter. The MOS expo- 
sures were taken in standard Full Frame mode, and for 
these files PATTERNS 1-12 were considered. The pn ob- 
servations were obtained in Extended Full Frame (EFF) 
mode. In view of the calibration uncertainties associated 
with the EFF mode, only PATTERN was considered, 
corresponding to the well-calibrated single events. 

The dedicated blank-sky event files accumulated by 
Lumb (2002) were used as background files. Events were 
extracted from each event file according to the same 
PATTERN selection as used for the source observations. 
The background files were then recast into the coordinates 
of the source observation, ensuring that source and back- 
ground products are extracted from the same regions of 
the detector, thus minmising detector variations. 

To ensure that all output products (spectra, surface 
brightness profiles, images) are corrected for vignetting, 
the photon weighting method of Arnaud et al. I|2001[l was 
used. This is applied to both source and background data 
sets. Note that while the background component induced 
by cosmic rays is not vignetted, source and background 
event files are treated in the same way, so that the correc- 
tion factor is the same and thus cancels. 

2.2. Flare removal 

Inspection of the MOS light curves in the source-free [10- 
12] keV band (Fig. shows that the background level 
during this observation is not constant. The MOS [10- 
12] keV light curves are very similar to the pn [12-14] keV 
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Fig. 1. The MOS1 light curve in the [10-12] keV band. All 
frames in the shaded area were excluded because of flare activ- 
ity. The 2a cutoff level (dotted line) , established after exclud- 
ing the shaded region, and the frames used for the "severe" 
analysis (dot-dash line), are also shown. Further details in the 
text. 

light curve. These high energy curves show that the overall 
level seems to be slowly increasing, with some low-level 
flaring activity appearing near the end of the observation. 
The fluctuations are not large in these bands, but could 
potentially have an effect on spectra in the outer regions 
of the cluster, where the emission begins to be dominated 
by the background. 

However, inspecting the light curves at lower energies, 
it is obvious that, for this observation at least, any cri- 
terion applied in the high energy band would miss some 
flares which appear to reach peak amplitude in lower en- 
ergy bands. 

It is prudent to exclude all frames near the end of 
the observation (the last ~ 5 ks), where flaring begins to 
become important. To test the effect of imposing differ- 
ent flare cutoff thresholds, three event files were then pro- 
duced, filtered at different levels of severity, increasingly 
diminishing the final exposure times. The background 
event lists were treated using exactly the same criteria. 

1. In the "2cr" event list, using the [10-12] keV light 
curve, frames were excluded which were outside a 2a 
threshold according to the method outlined in Pratt & 
Arnaud (2002). The final exposure times are 23/16 ks 
(MOS/pn). 

2. There are some observations where low-energy flares 
exist which do not appear at high energies. Light 
curves were thus extracted in several bands. For the 
MOS detectors, these bands were [0.3-1.4], [2.0-5.0], 
[5.0-8.0], [8.0-10.0] and [10.0-12.0] keV. For the pn, 
the [8.0-10.0] keV band was ignored because of strong 
instrumental lines, but a further light curve was ex- 
tracted from events with energies between [12.0-14.0] 
keV. The 2a filtering described above was then applied 
to each light curve, generating 5 GTI files. The obser- 



vation was then filtered using the intersection of all 
these GTI files for a final "intersect" event list. The 
resulting exposure times are 18/12 ks (MOS/pn). 
3. The first ~ 13 ks of the observation appears to be rel- 
atively quiescent, and has a mean count rate similar to 
that of the background files. A "severe" event list was 
generated by discarding all frames after this period, 
leaving 13/8 ks (MOS/pn) of useful data. 

After masking extraneous sources, surface brightness 
profiles were extracted from each event list in several dif- 
ferent bands. These profiles were background subtracted 
and compared. They are identical within their errors apart 
from differences in normalisation. To assess the effect of 
the filtering on spectral products, spectra were extracted 
in annular regions, centred on the peak of the X-ray emis- 
sion from the cluster (this is is discussed in more detail be- 
low). The annular spectra from each event file were fitted 
with a MEKAL model at the redshift of the cluster, ab- 
sorbed with the galactic absorption. MOS1, MOS2 and pn 
files were fitted both independently and simultaneously. 
The profiles from each event file were again consistent 
within their errors. 

It thus appears that the different flare filtering levels 
have little or no effect on the main output products, likely 
due to the relatively weak nature of the fluctuations. It 
is nevertheless important to establish this, as incomplete 
flare screening has been suggested as a source of discrep- 
ancy between XMM-Newton and Chandra observations of 
A1835 (Markevitch I2002j> . In an effort to strike a delicate 
balance between enthusiasm and caution, the "intersect" 
event lists are thus adopted for the remainder of this anal- 
ysis. 

2.3. Background subtraction 

For each of the output products (spectra and surface 
brightness profiles), an equivalent product is extracted 
from the corresponding blank-sky event list. As discussed 
extensively in Pratt et al. ptJOTJI and Arnaud et al. <|2002[l . 
while the blank-sky observations provide an adequate de- 
scription of the hard X-ray and instrumental components 
of the XMM-Newton background, they do not necessarily 
represent the soft X-ray component of the background, 
because this is variable across the sky (see Snowden et al. 
HWjl . The method described in Pratt et al. lj2TM|) and 
Arnaud et al. I|2002(l is used to correct for the difference 
of the background at low energy. This correction is applied 
to all spectra and surface brightness profiles discussed be- 
low. 

The cosmic ray background is variable at the ~ 10%- 
level, and so it is frequently necessary to normalise the 
background to the level of the observation. The count rate 
in the [10-12] keV and [12-14] keV (pn) bands, for MOS 
and pn respectively, is used in this analysis, with each 
camera treated separately. 
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Table 1 . Results of analytical fits to the gas surface brightness 
profile, errors are 90% confidence. 
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Fig. 2. An XMM/DSS overlay image. The X-ray contours are 
from an adaptively smoothed [0.3-2.0] keV non-background 
subtracted MOS+pn image, and are logarithmically spaced by 
0.15 dex with the lowest contour at 4.2 x 10~ 3 cts s _1 arcmin 2 . 

3. Gas density distribution 

3.1. Morphology 

An XMM/DSS overlay image is shown in Fig.|2 The clus- 
ter exhibits morphologically symmetric X-ray emission, 
centred exactly on the early type galaxy [WCB96] ACO 
1983 B, at z = 0.04405 (Wegner et al. \FMl\TW$\i . Note 
that this galaxy is not the Brightest Cluster Galaxy as 
defined by Postman & Lauer l|1995fl . which is located 12' 
further away in the North, well beyond the X-ray emis- 
sion, but the second ranked galaxy: [PL95] ACO 1983 G2. 

3.2. Surface brightness profile 

Tha gas density profile can be derived from the sur- 
face brightness profile, knowing the emissivity in the 
considered energy band. This is however more tricky 
for cool clusters than for hot clusters. For clusters like 
A1983 (kT ~ 2 keV), the exponential cut-off of the 
brcmsstrahlung emission lies at low energies (E ~ kT). 
Furthermore, line emission is prominent, specially the FeL 
blend (see Fig. 0} and radiative recombination radiation 
contributes significantly to the continuum. The emissivity, 
A, in any considered energy band thus depends sensitively 
on the abundance and temperature. As will be seen below 
(see Sect. l4.3l and Sect. 14.6)) radial variations of these quan- 
tities exist in the cluster centre. Thus, the corresponding A 
profile has to be taken into account in the derivation of the 
gas density profile from the surface brightness. The fact 
that the temperature and abundance profiles are deter- 
mined with a cruder spatial resolution than the Sx profile 
introduces a second complication and is a potential source 
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of systematic uncertainties. The surface brightness profile 
was thus extracted in the [0.3-3.0] keV band with the [0.9- 
1.2] keV band excluded to minimise the contribution from 
the FeL blend, the feature which is the most sensitive to 
abundance and temperature. 

For each camera, the azimuthally averaged surface 
brightness profile of the cluster and corresponding back- 
ground was produced, after masking of bright sources 
identified in the SAS pipeline detection. The profiles were 
made directly from the event files by binning weighted 
events into circular annuli about the emission peak. The 
background subtraction, performed as described above in 
Sect 12.31 was undertaken separately for each camera. The 
resulting profiles are in excellent agreement with each 
other (bar the expected normalisation factor), and so the 
MOS and pn profiles were coadded. The total profile was 
then rebinned so that a S/N ratio of 3cr was reached. 

The background subtracted surface brightness profile 
was then corrected for the remaining radial variation of 
A. For that purpose a parametric model (in the form of 
Eq. [5J was used for the derived abundance and temper- 
ature profiles. The corresponding emissivity profile A(#) 
was then computed in the considered energy band using 
an absorbed thermal model convolved with the instrumen- 
tal response. The observed surface brightness profile was 
then divided by A(8) normalised to its value at large radii. 
This correction is about 50% for the central bin, but drops 
to about 10% at ~ 0.3' ~ 20 kpc, and is negligable be- 
yond ~ 1'. 

This corrected Sx profile is shown in Fig. [3J Cluster 
emission is significantly detected out to 8'. 4 or ~ 600 kpc. 
It is now directly proportional to the emission measure 
profile, EM(r) = f dl, where n c is the electronic 
density, which can be derived using the emissivity at large 
radii. 



3.3. Gas density profile 

It is convenient, especially for the total mass determina- 
tion below (Sect.EJ), to have an analytical description of 
the gas density profile at all radii. The surface brightness 
profile S{8) was thus fitted with various parametric mod- 
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Fig. 3. The combined EPIC surface brightness profile of 
A1983 in the [0.3-3.0] keV band, with the [0.9-1.2] keV band 
excluded and corrected for the dependence of the emissivity, A, 
on the temperature and abundance. The profile is background 
subtracted and corrected for vignetting. Dot-dash line: best fit 
/3-model to the outer (6 > l'.O) cluster region. The black and 
red lines are the best fit BB and KBB models convolved with 
the XMM-Newton PSF. See Sect. for model details and 
Table for the best-fit parameter values. 



els, all of which were convolved with the XMM-Newton 
PSF and binned in the same way as the observed profile. 

A standard /3-model is a good description of the outer 
regions, but fails nearer the centre (8 < 1'), where the 
data show a slight excess above that predicted from the 
fit. Progressively cutting the central region decreases the 
reduced x 2 . For the single /3-model, the % 2 becomes stable 
after the inner ~ 1' is excised; the fit results give j3 = 
0.65 ± 0.04 and r c = l'.Sl^JJ, with X 2 = 78.52/60 d.o.f. 

Pratt & Arnaud ( 2002 ) discuss alternative parameter- 
isations for cases where a central excess is seen. These 
data were thus fitted with their double isothermal /3-model 
(BB) and the generalised double /3-model (KBB); the lat- 
ter allowing a more centrally peaked gas density profile 
The resulting best fit models are given in Tabled Note 
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Fig. 4. The global spectrum of A1983: MOS1 (black), MOS2 
(red) and pn (green), shown with the best fitting combined 
model at kT = 2.13 ± 0.08 keV with Z = 0.40 ± O.O5Z . 

that the profile obtained beyond i? cut for the BB and KBB 
models is a classic /3-model. 

The BB and KBB models can be compared using an F- 
test. The F-test indicates that the KBB model is a better 
fit at the 42% confidence level, suggesting that the density 
distribution in the core is not significantly more peaked 
than for a conventional /3-model. The BB model is thus 
adopted for the remainder of this analysis. 

4. Temperature and abundance distribution 

4.1. Preliminaries 

In the following, source and background spectra were 
extracted from the weighted event lists. In the case 
of the pn, these spectra were also corrected for out 
of time events, although in practice this makes little 
difference to the results. The spectra are fitted between 
0.3 and 6.0 keV (cluster emission is not detected above 
background beyond the upper limit), and the following re- 
sponse matrices were used: ml_thinlv9q20t5r6_all_15.rsp 
(MOS1), m2_thinlv9q20t5r6_all.15.rsp (MOS2) and 
epn_ef20_sY9_thin.rsp (pn). 

4.2. Global spectrum 

As a consistency check between the three cameras, a global 
spectrum was extracted from each event file using all 
events within the radius of detection (8'. 4, as established 
in Sect. 13. 3|) . Each spectrum was fitted with an absorbed 
MEKAL model with the redshift fixed at z = 0.0442, 
the mean redshift of the cluster. After checking that the 
fit results agreed when the Nh was left free, the absorp- 
tion was fixed at the galactic value of Nh = 2.04 x 10 20 
cm~ 2 (Dickey & Lockman 11990(1 : the free parameters of 
each fit were thus the temperature and the abundance 
(with respect to the solar abundances given by Anders 
k Grevessep89 ). The global spectral fits yield 2.25±g;?| 
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Table 2. Projected radial temperature profile results. Errors 
are la, as shown in Fig. 
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keV, 2.06l°;^ keV and 2.06±£^ keV for MOS1, MOS2 
and pn, respectively, values which are in excellent agree- 
ment. The combined MOS/pn fit yields kT = 2.13 ± 0.08 
keV and Z = 0.40 ± O.O5Z . 

The spectra are well fitted down to 0.3 keV. This is 
in contrast to a previous XMM-Newton observation of 
A1413 (Pratt & Arnaud I2UTJ2|| . where, below ~ 0.6/1.0 
keV (MOS/pn) a marked excess was seen. This could in- 
dicate that the excess in A1413 may well be a true soft 
excess, and not due to residual calibration problems. The 
global spectrum of A1983 is shown with the best fit model 
in Fig. |1 



4.3. Radial temperature profile 

A radial temperature profile was produced by excluding 
point sources and extracting spectra in circular annuli cen- 
tred on the peak of the X-ray emission. The widths of the 
rings were chosen so that a minimum 5cr detection in the 
[2-5] keV band was reached. This was possible for all but 
the final annulus. A minimum width of 30" was also im- 
posed, corresponding to the 90% encircled energy radius 
of the MOS PSF. The spectra were binned to 3<r above 
background level to allow the use of Gaussian statistics. 

The spectra were fitted with the absorbed MEKAL 
model described in the previous Section. In all cases the 
MOS normalisations were tied together, but the pn nor- 
malisation was allowed to vary. In the final two annuli, 
the abundance was frozen at the abundance found for the 
third-last annulus. Spectra from different cameras were fit- 
ted individually at first: these results are identical within 
their errors, and so a further simultaneous MOS+pn fit 
was undertaken. The results are detailed in Table [21 and 
shown in Fig. [S] 
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Fig. 5. The total XMM-EPIC radial temperature profile of 
A1983. The solid line is the best fit to the temperature profile 
with a function of the form given in Eqn. [5] Errors are la 
confidence level. 

Jones & Forman (1997), who find evidence for a weak CF 
The cooling time can be calculated using 



of ~ 6M, 



© yr 



t 



cool 



2.9 x 10 10 yr 



kT 



1 keV V 10 



"H 



(1) 



from Sarazin 1*19861 With the central density derived from 
the BB model fit, t coo i ~ 8.5 x 10 9 yr. This value suggests 
that cooling, if indeed present, should be rather weak. 

Starting from the central annulus and working out- 
ward, the spectra were fitted with an absorbed two tem- 
perature MEKAL model. Except for the central [0'- 0'.66] 
spectrum, the addition of a second spectral component 
either does not provide a significant (as defined by appli- 
cation of the F-test) reduction in the x 2 , or, more often, 
the second component parameters are completely uncon- 
strained. 

The [0'- 0'.66] spectrum was then fitted with vari- 
ous models as detailed in Table |21 The best-fitting model 
is (very marginally) the absorbed MEKAL+MKCFLOW. 
The same spectrum fitted with a single temperature plus 
cooling flow model with M fixed to 6Mq yr -1 (as ob- 
tained by White et al. I1997[) yields xl — 1-93 and shows 
significant residuals below 1 keV. 



4.4. Is there a cooling flow? 

The slight and gradual decrease in temperature towards 
the centre indicates that there might be a cooling flow 
(CF). Indeed this has already been suggested, on the basis 
of a deprojection analysis of the Einstein data by White, 



4.5. Projection and PSF effects 

PSF and projection effects will be a particular problem 
if radial bins smaller than the 90% encircled energy ra- 
dius of the PSF are used, or if the surface brightness of 
the cluster is very peaked (i.e., the luminosity is concen- 
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Table 3. Comparison of spectral fits to the [0'- 0'.66] spectrum. The Ah has been fixed to the galactic values as detailed in 
Sect 14.21 All errors are given at the 90% confidence level. 



Model 


kT 


Z 




M 






(keV) 


(Z Q ) 


(keV) 


(M e yr- 1 ) 


X 2 /d.of 


wabs(mekal) 
wabs(mekal+mekal) a 
wabs(mekal+mkcflow) b 


l-88_0 .14 
n i O + 0.19 
^•1»_0.16 
qq~(~0.32 


0.56±^ 
0.85«;i« 


0.82±g;^ 


q r, n + 0.80 


331.4/307 
330.4/305 
330.1/305 



a In this fit the abundance of heavy elements is linked between the two components 

b The maximum temperature of the CF is tied to the temperature of the thermal component, and the abundances are linked 
between the two components. The low temperature cutoff of the CF is found to be 0.63lo'.2o 

keV. 
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Fig. 6. The radial abundance profiles for Fe (black) and Si 
(grey). 



4.6. Abundance gradients 

There is strong line emission from several elements at the 
average temperature of A1983, however, significant radial 
constraints could only be put on Fe and Si abundances. 
(The Mg K line falls in the energy band affected by the 
Al K and Si K fLourescent energy lines in the MOS back- 
ground, and so was ignored in the present analysis.) 

The radial abundance profile of Fe shown in Fig. was 
obtained from an absorbed VMEKAL model fit with all 
a-processed elements and Fe treated as free parameters. 
To augment the statistics outside the centre, the spectra 
were reaccumulated in annuli twice the width as those 
used for the temperature determination. The Fe gradient 
is relatively flat except for the central annulus, where there 
is a pronounced jump from Q.2>Zpe,Q to 0.5Zp e ,Q- 

The spectra had to be reaccumulated in bins twice 
as large again for constraints to be found on the Si pro- 
file, shown in the same Figure. The Si abundance jumps 
from 0AZsi,o to 0.8Zs ijQ in the central annulus, but out- 
side this region the Si abundance is essentially constant. 
Because of the large errors on the Si abundance, it is dif- 
ficult to obtain useful constraints on the radial variation 
of the Si/Fe abundance ratio. 



trated towards the centre), or if there are significant tem- 
perature gradients. None of which is true for A1983: the 
radial temperature profile, extracted in bins with width 
greater than the 90% encircled energy radius of the PSF, 
is relatively flat out to the limit of detection, and, while 
the surface brightness profile does show a slight peak to- 
wards the centre, this is not a large effect. In this context 
A1983 can be compared directly with A1413, a cluster 
with a relatively flat temperature profile and a peaked sur- 
face brightness distribution. The fully PSF-corrected and 
deprojected temperature profiles of A1413 are consistent 
with the projected profile within the la errors. Since de- 
projection and/or PSF correction are likely to have little 
effect on the temperature profile, except to decrease the 
S/N, the projected profile is used throughout this analy- 
sis. 



4. 7. Hardness ratio image 

A hardness ratio (HR) image offers a relatively simple way 
of analysing the two-dimensional temperature structure of 
the cluster, and to this end, source and background images 
were produced in the [0.3-0.9] keV and [2.5-6.5] keV bands, 
chosen to avoid strong line emission. Exposure maps were 
generated for each image in each band using the SAS task 
eexpmap. 

To avoid large numbers of negative pixels, especially in 
the external regions, it is preferable to smooth images be- 
fore undertaking any mathematical operations. The SAS 
task asmooth allows a smoothing template to be defined, 
ensuring that each region of each image can be smoothed 
in exactly the same fashion. Since the high energy im- 
ages have a lower sigal to noise ratio, the template is best 
defined using these images. For the present analysis, a 
template was defined by adaptively smoothing the MOS1 
[2.5-6.5] keV image, and this template was then applied 
to each source and background image. 
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Fig. 7. The hardness ratio map of the cluster extending to 
a radius ~ 7' from the centre, from all EPIC cameras. The 
dynamic scale has been chosen to emphasise structure, is linear, 
and ranges from -0.75 (black) to -0.375 (white). The lines are 
due to the gaps between CCD chips. 

The corresponding background was normalised based 
on the effective exposure times and subtracted from each 
image. The background subtracted images were then di- 
vided by their respective exposure maps, and the hard- 
ness ratio image was calculated for each camera using HR 
= (image[2. 5-6.5 keV] -image [0.3- 1.4 keV])/ (image [2. 5- 6. 5 
keV]+image[0.3-1.4 keV]). The MOS and pn results are 
in excellent agreement, and the combined MOS/pn HR 
image is shown in Fig [7] 

In deriving these results, the difference between the 
blank-sky background and the A1983 background at low 
energy has been neglected. This would only affect the low- 
energy images as there is essentially no difference above 
~ 2 keV, and would preferentially affect the outer regions 
where the background becomes a significant fraction of the 
detected flux. In not taking this into account, the absolute 
HR values will be incorrect. However, the overall temper- 
ature structure will still be well described. Tests with the 
annular spectra (discussed above) show that the difference 
between the backgrounds begins to play a role at a radius 
of about 6', making the uncorrected spectra seem hotter 
than they actually are by about 10%. This would explain 
why the external regions of the map seem slightly harder. 
A direct translation of HR values into temperatures has 
not been attempted because of the complexity properly of 
treating this problem. 

The HR image shows quite a lot of structure, with 
point sources to the NE and SW of the centre standing 
out clearly. There appears to be a cool elongation run- 
ning approximately SE-NW, which has a radial extension 
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Fig. 8. The gas entropy (S = kT/n^ 3 ) profile derived from 
the temperature profile and the double /3-model (BB) for the 
gas density. The dotted line represents the entropy profile ob- 
tained using analytic models for the gas density and temper- 
ature distributions. The dot-dash line is the entropy obtained 
in the central bin if the spectrum is modelled with a two phase 
gas and the hotter of the two components is used to calculate 
the entropy. 

of ~ 2' (roughly 140 kpc). However, it should be stressed 
that the colour table has been adjusted to emphasise any 
structure, and in fact, the hardness ratio variations are 
not large at all within the region covered by the map (as 
can be seen from the colourbar). As a rough guide, the 
variations in hardness ratio present in this image repre- 
sent temperature variations of order 20%, a value which is 
entirely consistent with the radial values seen in Sect. 14.31 

5. Entropy profile 

The gas entropy was determined from the BB analyti- 
cal model fit to the gas density profile and the observed 
temperature profile, with entropy S = kT/n^ 3 . The tem- 
perature profile was also modelled with a function of the 
form: 

T = T +T 1 [(r/r c )V(l + (r/rc)")] (2) 

(c.f., Allen et al. ISUUTIl with T = 1.86 keV, T x = 0.41 
keV, r c — 1'.46 and r\ — 1.75, and the entropy calculated 
using this profile. The resulting entropy profiles are shown 
in Fig. |S| The uncertainty in the entropy profile is domi- 
nated by the uncertainty on the temperature distribution. 
Typical errors, corresponding to the error in each bin of 
the temperature profile are indicated in the figure. As dis- 
cussed above in Sect. 14.41 there is some evidence that the 
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Fig. 9. Integrated total gravitating mass and gas mass distri- 
butions. The total mass is shown with la errors. The dashed 
line shows the mass profile calculated from the analytical tem- 
perature and gas density distributions. The solid line through 
the points is the best fit NFW profile, with c = 3.75 and 
?"200 = 1480. In this and subsequent Figures in Section |S| this 
best-fitting NFW mass model is used to compute the scaled 
radius given on the top axis. 



central bin may have more than one temperature compo- 
nent. The entropy in this central bin calculated with the 
hotter of the two temperature components is also shown 
in the Figure, it agrees with the single temperature values 
within their errors. 



6. Mass analysis 

6.1. Total mass profile 

6.1.1. Calculation of the profile 

The total gravitating mass distribution shown in Fig. [§] 
was calculated under the usual assumptions of hydrostatic 
equilibrium and spherical symmetry using 



M(r) 



kT r 
G/im p 



d In i 



dim 



dim 



(3) 



where G and m p are the gravitational constant and pro- 
ton mass and /i = 0.609. The profile itself was calculated 
with an adapted version of the Monte Carlo method of 
Neumann & Bohringer (1995), which takes as input the 
analytical parameters for the gas density profile and the 
measured temperature profile. The original method calcu- 
lates random temperature distributions within the error 
bounds of the observed temperature profile, using a 'dif- 
fusive' process characterised by two parameters, a win- 
dow size and a step parameter. A significant advantage 
of the method is that it guarantees a certain smoothness 
of the temperature profiles and thus naturally limits 'un- 
physical' profiles with large oscillations. However, for this 
cluster, probably due to the fluctuations observed in the 
temperature profile (see also below), it was found difficult 
to choose objectively the window and step parameters, on 
which depends the final dispersion on the derived mass 
profiles and thus the error on the mass profile. To avoid 
this problem, a simpler, although less efficient, way to gen- 
erate random temperature profiles was used. A random 
temperature at each radius of the measured temperature 
profile was simply generated assuming a Gaussian distri- 
bution with sigma equal to the ler error, and a cublic spline 
interpolation was used to compute the derivative. As in 
Neumann & Bohringer (1995) only 'physical' temperature 
profiles are kept, i.e. those yielding to monotically increas- 
ing total mass profiles. In total 1000 such profiles were 
calculated. 

To include the error due to the uncertainty on the mod- 
elling of the gas density profile, the error on the density 
gradient (i.e., dlnn g /dlnr) was calculated. For this, the 
surface brightness profile was fitted at each radial point 
with dlnrig/dlnr as the free parameter, enabling the er- 
rors on this quantity to be calculated from the \ 2 varia- 
tion, the other parameters (normalisation, r c in , r c ) having 
been optimised. The final uncertainty values for the mass 
points are derived from the quadratic addition of these 
errors with those due to the temperature profile (output 
from the Monte Carlo method). The Monte Carlo mass 
profile, with appropriate uncertainities, is shown in Fig. El 

The mass profile can also be derived analytically using 
the temperature profile model as given in Eqn.|2]and the 
best fit BB model. This result, shown as the dashed line 
in Fig. is in excellent agreement with the results from 
the Monte Carlo analysis. 
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Fig. 10. The integrated gas mass fraction as a function of 
radius. Errors are calculated from the quadratic combination 
of la errors on the total mass and the gas mass. The solid line 
is the gas mass fraction derived from the gas mass and the best 
fitting NFW mass profile. 



6.1.2. Modelling of the mass profile 

The total mass profile was fitted with the Navarro, Frenk 
& White l(TM7|> profile, given by: 



Pc(z)S c 



{r/r s )(l+r/r s y 



(4) 



where p(r) is the mass density and p c {z) is the critical 
density at the observed redshift, which, for a matter dom- 
inated f2 = 1, A = Universe is: 



3H| 
8ttG 



(i + *)<*; 



(5) 



and where 5 C can be expressed in terms of the equivalent 
concentration parameter, c: 



Sr. = 



200 



3 [ln(l + c) - c/(l + c)} ' 



(6) 



The radius corresponding to a density contrast of 200 is 
^200 = cr s . 

This produced a fit with % 2 = 9.14/7 d.o.f, which is not 
a bad fit. The best fit NFW parameters are: c = 3.75 ± 
0.74, r s = 394 ± 124 kpc, giving r 200 = 1480 kpc and 
M 200 = 2.15 x 10 14 M . This best-fit NFW model is shown 
overplotted on the mass profile of the cluster in Fig. [9] 
Note that in this and other Figures in this Section, the 
virial radius (7200 ) shown is that calculated from the NFW 
fit. 
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Fig. 11. The integrated radial mass-to- light ratio. Errors are 
derived quadratically from the mass and luminosity errors in 
each bin. The solid line is the mass-to-light ratio determined 
from the best-fitting NFW mass model. 



The alternative Moore et al. (1999) profile was also 
tried, which is given by: 



p(r) 



where 



Pc{z)5 c 



(r/r,)3/a i + ( r /r s f 2 



100c 3 



lnfl 



„3/2\ 



(7) 



(8) 



With all parameters free, this profile gives a worse fit than 
the NFW (x 2 = 12.22/7 d.o.f), and moreover the best-fit 
scale radius in this case is a rather unlikely 55 Mpc. The fit 
was repeated with a upper limit to the scale radius fixed 
at twice the best-fit NFW value (i.e., 788 kpc). Here the 
fit is worse again (\ 2 — 17.09/7 d.o.f), but gives similar 
r2oo and A/200 values to the NFW fit. The best-fit scale 
radius value is at the maximum allowed by the fit. In the 
light of these results, it seems likely that these data do not 
have the required radial reach to put useful constraints on 
this model. 

6.2. Gas mass profile and gas mass fraction 

Figure [5] also shows the radial gas mass distribution de- 
rived from integration of the best-fit BB model. The vari- 
ation in the gas mass fraction with radius is shown in 
Fig. El This figure shows that the gass mass fraction rises 
rather abruptly in the first 200 kpc, beyond which it sta- 
bilises at a value of ~ 7%. 



6.3. Mass to light ratio 

The extensive optical observations of A1983 can be used to 
derive the radial mass-to-light ratio, which is an indication 
of the extent to which the light traces the mass. 
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The contribution of the Central galaxy (cG) was the 
first to be estimated. The photometric data of Saglia et 
al. 1|1997[) was used: a blue magnitude of B = 16.0 mea- 
sured in an aperture of 8 ap — 14.95", together with a 
surface brightness profile following a de Vaucouleurs law 
with an effective radius of 6 C — 14". A total B magnitude 
of B = 15.3 was derived, which corresponds to a total blue 
luminosity of Lb.cG = 8.8 x 10 10 h^~ 2 L for a galactic 
extinction of A B =0.112 (Schlegel et al. HTM)> . The lu- 
minosity within any radius of interest was estimated using 
the aproximation of Mellier & Mathez l|1987fl for depro- 
jecting the de Vaucouleur surface brightness profile. 

Girardi et al. (2002) derive a total Bj luminosity of 
A1983 from APS data. After correcting for the complete- 
ness limit they obtain L Bj = 5.97/7.44 x 10 11 L Q within 
0.99 Mpc for a Hubble constant of 100 km s" 1 Mpc" 1 . 
The two luminosity values are derived using different pro- 
cedures for the fore/background correction, so the average 
of the two values was used. The luminosity was then con- 
verted to the B band assuming Z^/Lb^o = 1.1 Lb /Lb,© 
(Girardi et al. I2002|) . and the luminosity of the cG sub- 
tracted. The resulting luminosity (excluding the cG) is 
then Lb = (2.4 ± 0.2) x 10 12 h^ 2 L Q within a projected 
radius of 1.98 tuj~ Mpc. A modified King profile with core 
radius of r c = 20 hijj, 1 kpc and a = 0.61 (Girardi et al. 
1995} can then be used to model the galaxy distribution 
and thus interpolate the luminosity to the radius of inter- 
est. The 68% errors on the profile parameters are large: 
r c < 100 h^ 1 kpc and 0.5 < a < 0.77 (M. Girardi, private 
communication), and are the main sources of uncertainty. 
This luminosity profile was then added to the cG luminos- 
ity profile to derive the total cluster luminosity profile. 

The final total mass-to-light ratio is shown in Fig. 1111 
where the errors are computed quadratically from the 
mass and luminosity errors in each bin. 



6.4. Iron mass over luminosity ratio 

The Iron Mass over Luminosity ratio (IMLR = M-p e /L-Q) 
is a fundamental quantity for the understanding of ICM 
enrichment (e.g. Renzini 1997). Using the Fe abundance 
measured at large scale ([Fe/H] ~ 0.3 with typically 
15% error) the IMLR within 0.3 r 200 is IMLR = (5.3 ± 
2.2) x 10~ 3 hgQ 172 M /L . An overall IMLR was also de- 
rived by extrapolating the gas mass to r2oo and assum- 
ing that the abundance remains constant. In this case, a 
great deal of extrapolation is needed, and caution is ad- 
vised in the interpretation of the results. On the other 
hand, the luminosity estimate within that radius is more 
robust, being less sensitive to the assumed light distri- 
bution. The overall IMLR obtained using this method is 
IMLR = (8.9 ± 1.8) x 10~ 3 h^ 172 M /L . 



7. Discussion 

7.1. Mass profile: comparison with theoretical 
expectations 

There is now converging evidence that the mass profile of 
relatively hot clusters follows the cusped distribution pre- 
dicted from numerical simulations of dark matter collapse 
(David et al. 127)1711 Arabadjis, Bautz & Garmire l2T)T)2l 
Allen et al. I27)0"T1 Pratt & Arnaud EES}- Our results sug- 
gest that this is also true for cool clusters. 

The NFW profile is a relatively good fit to the A1983 
data, but not as good as for instance for A1413 (Pratt 
& Arnaud 2002). However, there are a several points of 
note. For the calculation of the mass profile, the radial 
temperature and gas density distributions have been used. 
The first and most obvious point is that the temperature 
distribution is not smooth. The radial temperature data 
used for this mass analysis exhibit fluctuations on the or- 
der of ~ 20%, and this will inevitably have an effect on 
the quality of the NFW fit to a mass model derived us- 
ing these data. Secondly, there is a point of inflection in 
the mass profile around 170 hjjjj 1 kpc, corresponding to the 
boundary between the two [3 models used to fit the surface 
brightness profile. In the future, our intention is to develop 
more robust models for the surface brightness modelling, 
but that is beyond the scope of this paper. 

The theory of cluster formation predicts that the con- 
centration parameter should be higher for lower mass sys- 
tems (e.g. NFW). The concentration parameter of A1983, 
c = 3.75 ±0.74, is lower than would be expected for a clus- 
ter of this mass, c ~ 9 (e.g Eke, Navarro & Frenk©98). 
It is also significantly lower than the value c = 5.4 ± 0.2 
derived for A1413, a cluster about 5 times more massive 
(Pratt & Arnaud l2002[) . However, studies of large numbers 
of simulated haloes by Bullock et al. l|200l| and Thomas 
et al. 1)200 l|l has shown that there is a wide dispersion 
in the value of c, so the result is perhaps not surprising. 
In addition, the original NFW profile was derived from 
the mass distributions of equilibrium haloes; in fact NFW 
specifically considered earlier epochs for some of their sim- 
ulated haloes, assuring equilibrium by excluding the most 
recently accreted subunit. As shown by Jing (2000), in a 
study of the diversity of density profiles of a large num- 
ber of massive haloes, the quality of the NFW fit depends 
on whether the halo is in equilibrium. In particular, unre- 
laxed clusters yield concentration parameters lower than 
the NFW values. Observationally, a close look at the X-ray 
image of A1983 (Fig. [21 indicates that, while the isophotes 
are regular at large scale, there is an extension towards the 
NE in the central regions. Moreover, as the X-ray temper- 
ature analysis has shown fSect. l4~7jl . there are radial and 
non-radial temperature variations of ~ 20%. Lastly, pre- 
vious optical observations (Escalera et al. 1994; Girardi 
et al. 11997(1 have found evidence for several substruc- 
tures in the field of this cluster, at least one of which 
appears to cover the central region. These observational 
points seem to suggest that A1983 is not as relaxed as ini- 
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tially it may appear, presumably due to an earlier merging 
event which may still be detectable in the X-ray tempera- 
ture and galaxy distributions. This shows the importance 
of knowing in detail the dynamical state of a cluster for 
the interpretation of the mass profile. 

Large simulations, like the Hubble Volume simula- 
tion, indicate a tight correlation between M200 and the 
dark matter velocity dispersion in the form of a = 
1075 [(E(z) M 200 /(10 15 h^ Mo)] 1 / 3 km s" 1 (Evrard & 
GioiaESEll, where E{z) = H(z)/H . The best fit NFW 
mass model, M200 = 2.15 x 10 14 h" 1 M Q , gives a ~ 520 
km s _1 . This is consistent with the robust, optically- 
derived value of the galaxy velocity dispersion a — 551+4 7 
km s _1 . This excellent agreement indicates: i) that the to- 
tal mass estimate is probably not too much in error and 
that the cluster cannot be too far from equilibrium, and 
ii) gives further support to the current modelling of the 
dark matter collapse. 

Finally, the normalisation of the M200-T relation for 
this cluster, defined as M 10 (200) = E(z) M 200 T w 3/2 , 
where Tiq is the temperature in units of 10 keV, is 
Mio(200) - 2.3 x 10 15 h^ 1 M Q . This is about 25% lower 
than the normalisation of Evrard et al. I|1996[) . We thus 
confirm, at a lower mass, the offset between the observed 
and theoretical normalisation of the M-T relation, al- 
ready observed for hot clusters (Allen et al. 120011 Pratt 
& Arnaud Enn3) ■ 

7.2. Mass versus light as a function of scale 

At large radii, the light seems reasonably to trace the 
mass. Within the relatively large errors, the mass-to-light 
ratio flattens beyond ~ 0.2 r2oo, reaching a value of 
M/L B = (135 ± 45) h 50 M /L Q at 0.3 r 200 (Fig. EJ. 
The value extrapolated at r 2 oo, using the best fit NFW 
mass profile, is similar: M/Lb = 150 I150 M0/L0, for a 
luminosity within that radius of Lb ~ 1-4 x 10 12 h^~ 2 L 
(estimated as described in Sect. 16.311 . 

This mass-to-light ratio value for a cluster of such a 
low luminosity is consistent with the general trend ob- 
served by Girardi et al. l|2002|l from poor groups to rich 
systems. Using a very large sample, these authors found 
a significant increase of M / L with luminosity with a 
best fit power law of M/M = 10- 1596 (i B /L ) 1338 
for Ho = 100 km s _1 Mpc -1 . Rescaling their results for 
I150 = 1 , and using the luminosity estimated at f2oo j leads 
to a mass-to-light ratio of M/L B ~ (100±20) h 50 M /L , 
taking into account the typical dispersion around the cor- 
relation. This is within the error bars of the observed 
value, although somewhat on the low side. 

On the other hand, David, Jones & Forman (|l995), 
who used X-ray mass estimates for a much smaller sample 
of 7 groups and clusters, did not find any evidence of mass- 
to light variation with mass. All objects in their sample fall 
in the range M/Ly = 100 — 150 hso M /L . Similar val- 
ues were derived by Cirimele, Nesci & Trevese ©97 ) from 
a sample of 12 clusters in the temperature range from 2 to 



10 keV. Assuming (Ly/Lv,©) = 1-3 (Lb/Lb,©), relevant 
for early type galaxies, this translates into a mass-to-light 
ratio in the B band of M/L B = 130 - 200 h 50 M /L . 
It must be noted that the mass-to-light ratio of A1983 
is also consistent with the David et al. 1)1995(1 finding. 
Clearly larger homogeneous X-ray samples, of size simi- 
lar to those now available in optical, are needed to see if 
there is a real contradiction between M/L trends derived 
respectively from X-ray and optical mass determination 
methods. 

The mass-to-light ratio decreases towards the central 
part of the cluster, where the light is dominated by the 
central galaxy fFig. lll|) . The mass-to- light ratio in the cen- 
tral bin, located at 23 kpc or about 1.4 times the effective 
radius of the central galaxy (Sect. EU, is M/L-q = 9 ± 2. 
This value is consistent with the typical mass-to-light ra- 
tio of the stellar population (e.g. Gerhard et al. 120011 and 
references therein). This suggests that there is little if any 
dark matter component in the central part of the galaxy. 
At r ~ 2.8r e the M/L-q ratio has increased to a value of 
~ 30. Interestingly, similar values and trends have been 
observed both around M87, the central very bright galaxy 
in the Virgo cluster (Matsushita et al. 2002J) and in iso- 
lated elliptical galaxies (Gerhard et al. I2001f) . two very 
different environments. 

7.3. Scaling law tests 

7.3.1. Testing the M gas -T relation 

The surface brightness profile of A1983 can be used to 
test the M gas — T relation. The surface brightness profile of 
the cluster was converted into an emission measure (EM) 
profile using 

£ M M = ^yW ; r = ^ } ,0 (9) 

where A(T, z) is the emissivity, taking into account the 
interstellar absorption and the spectral response, and r = 
cIa(z) is the angular distance at redshift z. 

This EM profile was then scaled according to the self- 
similar model, using the standard scaling relations of clus- 
ter properties with redshift and temperature 1 . In the first 
instance, the profile was scaled using the standard relation 
M gas cx T 15 . In the left-hand panel of Fig.^J this scaled 
EM profile is compared to the mean scaled EM profile of 
the ensemble of hot (kT > 3.5 keV), nearby (z < 0.06) 
clusters observed by ROSAT and considered by Neumann 
& Arnaud l|2001|) . It can clearly be seen that the form of 
the profile is similar to that observed for hotter clusters. 
This is interesting, as is implies that the gas distribution 

1 In all the Figures in this Section, and for easy comparison 
with previous works, the virial radius (r^oo), is calculated using 
the normalisation given in Evrard, Metzler & Navarro 1996 
viz., r 20 o = 3.69 [T/(10 keV)] 1/2 (1 + z)" 3/2 Mpc. Note that 
when comparing clusters of different temperatures, the exact 
value of this normalisation does not matter. 
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Fig. 12. The scaled emission measure (EM) profile of A1983 (data points) compared to the mean scaled EM profile of the 
ensemble of hot clusters (kT > 3.5 keV) observed by ROSAT and considered by Neumann & Arnaud 12UUH . The latter 
is indicated by the dashed line; dotted lines show the dispersion. On the left, the EM profile of A1983 is scaled using the 
self-similar M gas oc T 1 ' 5 scaling. On the right, it is scaled using the empirical M gas oc y 194 scaling of Neumann & Arnaud (2001 ) 



of this poor cluster is not more inflated than that ob- 
served for hotter clusters. In the BB model (Sect. ET31l the 
outer regions are modelled using a standard /3-model: the 
best-fit (3 value is f3 = 0.76±g;J|, a value which is roughly 
consistent with (but maybe a little higher than) the aver- 
age value for hotter clusters, (3 — 0.67, found by Neumann 
& Arnaud PUUT)) . 

However, there is an obvious problem with the normal- 
isation (there is a factor of two difference between the ref- 
erence profile and that of A1983). In the right-hand panel 
of Fig. El the EM profile has been scaled using the re- 
lation M gas oc T 1 - 94 . This relation, empirically-determined 
by Neumann & Arnaud H2001|) . was found significantly to 
reduce the scatter in the scaled profile of nearby hot clus- 
ters. It is clear that using this M gas —T relation, in other 
words, allowing a variable / ga s, has considerably improved 
the normalisation, allowing A1983 to fall within the dis- 
persion observed in the scaled EM profiles for hotter clus- 
ters. 

In summary, the gas distribution of this poor cluster is 
not more inflated than that observed for hotter clusters, 
and the scaling of the EM profile of A1983 is consistent 
with the steepening of the M gas —T relation observed for 
hotter clusters. 

7.3.2. The entropy profile and the S — T relation 

Here the XMM-Newton entropy profile of A1983 is com- 
pared with the XMM-Newton entropy profile of the hot 
(kT = 6.9 keV) cluster A1413 (for more details see Pratt 
& Arnaud 2 002fl . In a self-similar case, where / gas is con- 
stant, S oc (1 + z)~ 2 T . This scaling is shown in the the 
left-hand panel of Fig. El It can immediately be seen that 



the profiles are practically identical in form. This is itself 
is a remarkable result which will be discussed in detail 
later on. 

Ponman et al. (2003) have shown, with a sample of 66 
virialised systems, that the entropy measured at 0.1r2oo m 
fact scales as S oc 7 1 ~ - 65 . Implicit in this scaling is a de- 
pendence of the gas density n G , with system temperature. 
Applying this scaling to the entropy profiles of A1983 and 
A1413, shown in the right hand panel of Fig. El shows 
that the normalisation between the two profiles is con- 
siderably improved. This Figure also shows the S oc r 
behaviour, expected from analytical modelling of shock 
heating in spherical collapse (Tozzi & Norman EDDTJ and 
normalised to the scaled entropy of a 10 keV cluster from 
Ponman et al. (2003). It can be seen that overall, and 
down to surprisingly small radius, both profiles agree re- 
markably well with this prediction. It is also clear that, al- 
though the profiles flatten somewhat below about 0.1 r 2 oo, 
there is no evidence of an isentropic core. 

7.3.3. Discussion 

The role of non-gravitational processes and the breaking 
of self-similarity in clusters has long been a hot topic. 
Approaches have included analytic, semi-analytic and nu- 
merical modelling, including pre-hcating, internal heating 
by AGNs and SNs and/or cooling (see Voit et al. 120021 for 
a full review). 

A dependence of / gas on temperature has been pre- 
dicted in some numerical simulations (Dave et al. [5002), 
and, as Lx oc n\, has been suggested as a principal cause 
of the similarity break in the Lx~T relation. However, 
this variation has only been hinted at in the observations 
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Fig. 13. The scaled XMM-Newton entropy profiles of A1983 and A1413. The points show the data. The solid line running 
through the A1983 points is the entropy obtained when the temperature profile is modelled using Eqn.|2] the dotted line running 
through the A1413 points is the entropy obtained when the temperature profile is modelled by a polytropic model with 7 = 1.07 
(see Pratt & Arnaud 2002;for details). In the left hand panel, the profiles are scaled with the self-similar S oc (1 + z)~ 2 T scaling. 
In the right panel, the profiles are scaled using S oc (1 + z) -2 T ' 65 . The dashed line in the right-hand panel shows the S oc r 
behaviour expected from shock heating, normalised to the scaled entropy of a 10 keV cluster from Ponman et al. ( 2003 1 



(Mohr et al. 115551 : Arnaud & Evrard 1999 ). Both of these 
papers determined global / gas values using either /3-model 
or virial estimates for the total mass of the system. This 
approach, as stressed by the authors, is obviously model- 
dependent. The new XMM-Newton data allow the deriva- 
tion of much more precise / gas values, and allow the trac- 
ing of the / gas with radius. 

A1983 is obviously gas poor compared to the hotter 
cluster A1413. This can be seen not only in the raw / gas 
values (at 0.3r 200 , for A1983 / gas = 0.08 for A1983, while 
for A1413, /gas = 0.16) but also in the scaling behaviour 
of the emission measure profiles. Interestingly, Mushotzky 
et al. I|2003ll also find remarkably low / gas values in their 
analysis of XMM-Newton data from two galaxy groups 
(NGC2563 and NGC4325; gas mass fraction of - 4 - 6% 
at 0.3r2oo)- Thus there is strong evidence for some process, 
proportionally more important in low mass systems, which 
either biases intergalactic gas accretion, and/or removes 
gas from the hot X-ray emitting phase of the ICM after 
collapse. Note that the similarity of the profiles over a wide 
temperature range indicates that the observed variation of 
/gas is not simply due to a displacement of gas from the 
central regions to the outer regions. 

Closely related to the question of a variation of / ga s; 
and independent confirmation of that fact, is the entropy 
scaling behaviour. The relation S oc T°- 65 j observed by 
Ponman et al. 1(2003) acts to considerably improve the 
normalisation between the entropy profiles. This sort of 



modified S—T relation has been seen in numerical simula- 
tions incorporating radiative cooling Dave et al. ( 2002 ) . 

Perhaps the most interesting aspect of the A1983 and 
A1413 entropy profiles is their obvious similarity. There 
is no isentropic core and both profiles converge rapidly to 
the S(r) oc r value favoured by models of the shock- 
dominated regime. The scaled core sizes are exactly the 
same. 

As seen in e.g., Tozzi & Norman (|2001|l and Babul et 
al. ( 2002), preheating models suggest that in the low- 
est mass systems the accretion of the gas onto a potential 
well is isentropic, with the shock heating regime becoming 
dominant at higher masses. As a result the entropy profiles 
of the lowest mass systems have large isentropic cores, and 
these cores become progressively smaller with increasing 
mass because of the increasing importance of shock heat- 
ing. This is clearly inconsistent with the present results 
(as also seen in Ponman et al. 120031 and Mushotzky et al. 
I21TO . 

Given the evidence that preheating models cannot de- 
scribe the observed entropy profiles, one is forced to look 
towards models for which the observational tests are much 
less evident. These models invoke radiative cooling or 
heating by supernovae or AGN, or a combination of the 
two, as a mechanism for modifying the entropy of the ICM 
in lower mass systems. 

Pure cooling models, while yielding results which are 
qualitatively similar to the observations, seem to have 
problems with overcooling, which furthermore is affected 
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by the finite resolution of the simulations (Muanwong et 
al. 120021 Dave et al. l2002f) . It thus seems unlikely that pure 
cooling models can work without the addition of some 
form of feedback of energy into the baryonic component. 
This feedback may be due to star formation, or AGN heat- 
ing, or any other plausible but undiscovered mechanism. 
Moreover, there is a potential problem with the fate of 
the cooled gas, which must not exceed the observed mass 
fraction of the stellar component (e.g. Balogh et al. 120011 
Dave et al. E002I Borgani et al. l2T)0l|l . 

An interesting analytical treatment is that of Voit et 
al. f2002J), where the problem is attacked using a phe- 
nomenonological treatment of the entropy distribution of 
the gas. The ICM entropy is either truncated such that 
all the gas with a cooling time less than a Hubble time 
is removed (analagous to radiative cooling), or shifted by 
raising the entropy level throughout (analagous to pre- 
heating). A generalised radiative loss model is also used, 
which corresponds to a variable truncation of the entropy 
distribution (analagous to some form of cooling with feed- 
back). In common with other treatments, the preheated 
entropy distribution tends to show a larger isentropic core, 
and also the pure radiative cooling model seems to suffer 
from overcooling. The authors suggest that feedback is 
necessary to prevent this. 

It thus seems clear that for cooling to be a viable alter- 
native to the more popular preheating scenario, some form 
of heating/feedback is needed to prevent the overcooling. 
The question now is how to identify and distinguish be- 
tween the competing processes, but from an observational 
point of view it is unclear how this is to be done without 
more stringent constraints from modelling and numerical 
simulation. 

7.4. ICM enrichment 

It is instructive to compare the IMLR of A1983 with the 
ratio observed in hotter clusters. The results must however 
be treated with some care, in view of the heterogenity of 
the various analyses. 

From a compilation of the literature, Renzini (1997) 
found a relatively small scatter in the IMLR down 
to a temperature of 2.2 keV, with a average value of 
0.02 h^ 1/2 Mq/Lq. The IMLR then drops suddenly with 

— 1/2 

temperature, to IMLR < 0.01 h 50 ' M /L Q , for groups 
(kT < 1.8 keV). With an IMLR of 8.9 ± 1.8 10~ 3 (at r 200 ) 
for a temperature of 2.1 keV, A1983 lies at the boundary 
of the two regimes. It must be noted however the IMLR 
values considered by Renzini (1997) are not consistently 
estimated within the cluster virial radii. This is a poten- 
tial source of bias, since the IMLR is expected to increase 
with radius: whereas the Fe abundance profiles are flat at 
large scale (De Grandi & Molendi 12001), there is evidence 
that the gas distribution is more inflated than the galaxy 
distribution (e.g. Cirimele et al. H997fl . For instance, above 
3.5 keV, the data rely on the work of Arnaud et al. (jl992). 
These authors estimated the IMLR within a fixed radius 



of 3 Mpc. This corresponds to r2oo I0r a ~ 6.6 keV clus- 
ter. The IMLR ratios of cooler (hotter) clusters are thus 
estimated at a larger (smaller) fraction of the virial ra- 
dius and thus are probably over (under) estimated. It is 
beyond the scope of this paper to correct for this effect. It 
should however be emphasized that once this correction is 
made, one could well find a systematic increase of IMLR 
with temperature, across the whole temperature range, 
and that the IMLR of A1983 is clearly about two times 
less than for 'hot' clusters. 

An increase of IMLR from groups to rich clusters was 
observed by Finoguenov, Arnaud & David l|20011 Figure 8) 
for an IMLR estimated consistently at a fixed fraction of 
the virial radius (0.2 and 0.4 r 2 oo)- Unfortunately the large 
error bars and considerable scatter did not allow them to 
establish whether the IMLR systematically increases with 
temperature or tends to a constant above some tempera- 
ture limit. Again the IMLR ratio of A1983, 5.3 ± 2.2 10" 3 
at 0.3 T200; li es we ll within the observed trend and is again 
about two times smaller that the IMLR above 5 keV. 

On the other hand the abundance of 0.3 is similar 
to the abundance found in hot clusters (De Grandi & 
Molendi 2001). If confirmed, such an increase of IMLR 
with system temperature (i.e mass) while the abundance 
remains constant, will provide interesting information on 
the evolutionary history of the gas. 

A dependence on system mass of the Fe mass pro- 
duced by unit of stellar light is unlikely. The abundance 
is set by the ratio of the IMLR to the M gas /L ratio. It 
would require an unlikely conspiracy to keep the abun- 
dance constant: a similar increase of the IMLR and the 
M gas /L ratio, quantities which are obviously not driven 
by the same physics. Furthermore, such a variable IMLR 
is unlikely on theoretical grounds. In view of the abscence 
of abundance evolution with z, it is commonly assumed 
that the bulk of the Fe now seen in the gaseous phase 
was produced in galaxies at early epochs (probably before 
z ~ 2), and ejected through early galactic winds. Neither 
the heavy element production, determined by the chemi- 
cal evolution of the galaxies, nor the ejection efficiency are 
expected to vary with system mass. 

A more natural explanation of the observed variation 
of the IMLR is that we do not see all the Fe mass ever 
produced and ejected by the galaxies. This issue is actually 
closely related to the question of the variation of / gas with 
system mass. If the heavy elements ejected by galaxies are 
homogeneously distributed in the intergalactic medium, 
the X-ray determined IMLRx will scale with the fraction 
of the gas with is now in the hot phase of the virialised 
part of the IGM, i.e with the gas mass fraction, / gas , as 
determined from X-ray observations. On the other hand 
the abundance will rest constant. Whatever mechanism 
explains an increase of / gas with T, the increase naturally 
implies a corresponding increase in IMLRx- 

8. Conclusions 

The main conclusions are summarised below. 
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Results have been presented from the XMM-Newton 
observation of the poor cluster A1983 (kT = 2.1 keV), 
at z = 0.044. 

The gas density profile has been measured out to ~ 
8'. 4 and the temperature profile out to ~ 7', or ~ 
500 hjjjj kpc, corresponding to ~ 0.35 r2oo- The mass 
profile has been calculated out to the same distance 
assuming HE and spherical symmetry. 
The gas density profile is adequately modelled with a 
double isothermal /3-model (BB), in which the gas den- 
sity profile inside and outside a radius r cut is modelled 
by two different /3-models. The outer region /3-model 
parameter is (3 — 0.74. 

The temperature profile exhibits a small drop towards 
the centre, but remains roughly constant out to the de- 
tection limit. Both the radial and hardness ratio tem- 
perature information suggest fluctuations in temper- 
ature of < 20%. There is no evidence for multiphase 
gas except in the very central annulus, where a two- 
temperature or a cooling flow model gives a marginally 
better description of the data. 

The mass profile is consistent with an NFW profile 
with scale radius r s = 394 kpc and a concentration pa- 
rameter c = 3.75 ± 0.74. A Moore et al. I|1999[> profile 
is unconstrained, due to the lack of data at large radii. 
The low value of c may be indicative of the dispersion 
at lower masses observed in numerical simulations, or 
may suggest that the cluster is not completely relaxed. 
However, best fit NFW mass model gives a dark mat- 
ter velocity dispersion of a ~ 520 km s , in excellent 
agreement with the optically derived galaxy velocity 
dispersion of a = 551^4^ km s -1 . Finally, the normali- 
sation of the Af2oo — T relation for this cluster is ~ 25% 
too low with respect to the numerical simulations of 
Evrard et al. fflfflft. 

The gas mass fraction at 0.35r 2 oo is / gas = 8%. 
When scaled using the self-similar relation M gas oc T 15 , 
the emission measure (EM) profile of A1983 lies a 
factor of ~ 2 lower than the mean scaled profile of 
hot, nearby clusters discussed in Neumann & Arnaud 
(2001}, but has a similar shape. If the empirically de- 
termined relation, M gas oc T 194 , is used instead, the 
scaled EM profile of A1983 matches well the reference 
profile for hotter clusters. This is direct evidence for a 
variation of / gas with temperature. 
A direct comparison of the entropy profile of A1983 
with the XMM-Newton entropy profile of A1413 (kT = 
6.9 keV; Pratt & Arnaud 2002) shows that the scal- 
ing between the profiles is much improved when the 
empirically-determined relation S oc T is used, 
again supporting the suggestion of a variation of / gas 
with temperature. 

The entropy profiles of A1983 and A1413 are remark- 
ably similar in form, with no evidence of an isentropic 
core, suggesting that simple preheating models are not 
predicting the correct behaviour. This is in agreement 
with recent XMM-Newton and ASCA/ROSAT results, 



favoring models which change the ICM entropy by in- 
ternal heating, cooling, or a combination of the two. 

— The mass-to-light ratio rises from M/Lb = 9 ± 
2I150 Mq/L q at 23/i^q 1 kpc (typical of the mass-to- 
light ratio of the stellar population), to M/Lb = 
135 ± 451i5o M Q /L Q at 500htf kpc (0.35r 20 o; consis- 
tent with the trend of M/Lb with system mass noted 
by previous optical observations). 

— The iron mass over luminosity within 0.35 r2oo is IMLR 
= (5.3 ±2.2) x 10~ 3 h^ 1/2 Mq/L©. Extrapolating the 
gas mass to r 2 oo and assuming that the abundance 
remains constant, the overall IMLR is = (8.9 ±1.8) x 
10~ 3 h^ 1/2 Mq/L q . The IMLR is about two times 
less than the IMLR observed for clusters above 5 keV. 
This may also be connected to the variation of / gas 
with system mass. 

The results presented here have given some indication 
that the cluster population remains self-similar in shape 
down to low temperature, that the fundamental departure 
from the simplest self-similar model is an increase of / gas 
with temperature, and that there may indeed be some dis- 
persion in NFW fit parameters at the low mass end. Often 
in astronomy however, the prototype is the exception. In a 
forthcoming paper the results presented here will be tested 
with a larger sample of observations of poor clusters. 
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